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ABSTRACT 



^2 ' Aims. We report the discovery of a new white dwarf/M dwarf binary, HS 1857+5144, identified in the Hamburg Quasar Survey 
(HQS). 

Methods. Time-resolved optical spectroscopy and photometry were carried out to determine the properties of this new cataclysmic 
variable progenitor (pre-CV). 

Results. The light curves of HS 1857+5144 display a sinusoidal variation with a period of P or b = 383.52 min and peak-to-peak 
amplitudes of 0.7mag and 1.1 mag in the B-band and _R-band, respectively. The large amplitude of the brightness variation 
results from a reflection effect on the heated inner hemisphere of the companion star, suggesting a very high temperature 
of the white dwarf. Our radial velocity study confirms the photometric period as the orbital period of the system. A model 
atmosphere fit to the spectrum of the white dwarf obtained at minimum light provides limits to its mass and temperature of 
M wd ~ 0.6 - 1.0 M Q and T wd ~ 70 000 - 100 000 K, respectively. The detection of Hell A4686 absorption classifies the primary 
star of HS 1857+5144 as a DAO white dwarf. Combining the results from our spectroscopy and photometry, we estimate the 
mass of the companion star and the binary inclination to be M sec — 0.15 — 0.30 Mq and i ~ 45° — 55°, respectively. 
Conclusions. We classify HS 1857+5144 as one of the youngest pre-CV known to date. The cooling age of the white dwarf suggests 
that the present system has just emerged from a common envelope phase ~ 10 5 yr ago. HS 1857+5144 will start mass transfer 
within or below the 2-3 h period gap. 
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1. Introduction 

Post-common envelope binaries (PCEBs), i.e. detached 
white dwarf-main sequence binaries, originate from wide 
binaries comprising unequal main sequence components. 
Once the more massive star evolves through the giant 
phase and fills its Roche lobe, unstable mass transfer onto 
the unevolved star starts. The high accretion rate subse- 
quently brings the system into a common envelope phase 
(CE). Friction between the stellar components and the 
envelope shrinks the binary orbital separation and even- 
tually ejects the envelope from the system, resulting in the 
majority of PCEBs ha ving orbital periods of a few days 
(jWillems k Kolbll2004h . The binary separations of PCEBs 
are believed to be further reduced through angular mo- 
mentum loss via magnetic braking and/or gravitational 
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radiation. The main-sequence star will ultimately fill its 
Roche lobe and start mass transfer onto the white dwarf, 
which turns the binary into a cataclysmic variable (CV). 
While this general scenario is widely accepted, the details 
of the evolution through the common envelope phase, as 
well as the subsequent orbital angular momentum loss are 
poorly understood. Consequently predictions made by bi- 
nary population models are rather uncertain. 

As PCEBs are simple objects in terms of their stellar 
components, they offer a large potential in constraining 
and calibrating the physics of both CE evolution and or- 
bital angular momentum loss. This has implications for 
a wide range of astronomical objects such as SNIa pro- 
genitors, X-ray binaries, or neutron star binaries as likely 
progenitors of short gamma-ray bursts. However, until 
recently PCEBs received little observational attention, 
largel y due to the lack of a dedicat ed search for such sys- 
tems. ISchreiber k Giinsickel ( 2003 ) analysed the proper- 
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Fig. 1. A 10' x 10' finding chart for HS 1857+5144 ob- 
tained from the Digitized Sky Survey. The J2000 co- 
ordinates of the star are a = 18 h 58 m 32.12 s and 6 = 
+51°48'57.74". The comparison and check stars used in 
the photometry are marked by 'C1'-'C4' (see Sect. 12. 21 for 
details). 

ties of 30 well-studied PCEBs and showed that the known 
population of these systems is extremely biased towards 
young systems consisting of hot white dwarfs and late type 
companions, as the majority of the known PCEBs were 
initially selected as blue objects. A substantial improve- 
ment in the statistics of PCEB properties will be possi- 
ble through the exploitation of large extragal actic surveys 
such as the Hamburg Quasar Survey fHQS. [Hagen et aT 



1995) or the Sloan Digital Sky Survey (SDSS. lYork et al 
2000). 



Here we present a detailed follow-up study of the new 
PCEB HS 1857+5144, which has been discovered in our 
ongo ing effort to identify CVs and pre-CVs in the H QS 
(see Gansicke et al. l2002t Aungwerojwit et al. 2005 for 
details on the project). In Sect. [5] we describe the ob- 
servations and data reduction. The orbital ephemeris of 
HS 1857+5144 is determined in Sect.[3J In Sects.gHS] we 
analyse the nature of the stellar components and discuss 
the future evolution of the system. 

2. Observations and Data Reduction 

2.1. Spectroscopy 

An identification spectrum of HS 1857+5144 was ob- 
tained in August 1992 with the Boiler & Chivens 
Cassegrain spectrograph on the 2.2-m telescope at Calar 
Alto Observatory. The spectrum is characterised by a 
blue continuum superimposed by strong Balmer emis- 
sion lines. Subsequent time-series, intermediate resolu- 
tion spectroscopy of HS 1857+5144 was performed in 
July 2004 at the 2.7-m Harlan J. Smith telescope 
at McDonald Observatory equipped with the Large 
Cassegrain Spectrograph (LCS), covering ~ 10 h (35 



spectra) in total. The individual spectra were obtained 
through a l"slit and grating #43 and imaged on the 
800 x 800 pixel Til CCD camera. This setup provided 
access to the AA3670 — 5050 wavelength range at 3.5 
A spectral resolution. The reduction and optimal extrac- 
tion of the spectra were performed using standard long-slit 
spectroscopy packages within IRAF0. Four additional high 
resolution spectra were obtained using the Intermediate 
Dispersion Spectrograph and Imaging System (ISIS) on 
the 4.2-m William Herschel Telescope (WHT) on La 
Palma in July 2006, covering the orbital minimum, max- 
imum, and the quadrature phases. The blue arm of the 
spectrograph was equipped with the R1200B grating, pro- 
viding a spectral resolution of ~ 1 A over the wavelength 
range 4200 - 5000 A. The data were reduced using the 
Figaro within the Starlink package as well as Pamela 
and Molly written by T. Marsha 



Spectroscopic characteristics. Figures^] and [3] illustrate 
the orbital phase-dependent variation of the emission 
line profiles of HS 1857+5144 from the McDonald and 
WHT spectra, respectively. The emission lines greatly 
vary in strength over the orbital cycle, with maxi- 
mum line fluxes occurring at tp ~ 0.5 and disappear- 
ing around ip ~ 0.0. The dominant lines detected in 
the McDonald spectra belong to the Balmer series, but 
the high quality of the WHT spectra reveals a multi- 
tude of narrow emission lines, e.g. the Ciii/Nm A4650 
Bowen blend, Mgn, Nil, Oii/Ciii blend, and Niii/Fem 
blend. This emission spectrum is a characteristic of 
known PCEBs containing a cool secondary star ir- 
radiated by a very hot primary component, such a s 
BEUMa dFerguson et alJll98li iFerguson fc JamesHl9941 ). 
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(jSing et al.l 120041 ). The strengths of the Balmer emission 
lines vary in phase with those of Hei, Ciii/Nm, indicat- 
ing that all emission lines come from the same source. 
The WHT spectra also reveal that the Balmer emis- 
sion lines have double-peaked profiles, with peak separa- 
tions of ~ 170 km s -1 for H/3 (Fig. [31 small window) and 
~ 150 km s -1 for H7 which are most likely caused by non- 
LTE effects in t he strongly irradiate d atmosphere of the 
companion star ((Barman et alJl2004 ). 

During the orbital faint phase, ip ~ 0.0, weak absorp- 
tion Balmer lines are observed in both sets of spectra, with 
He 11 A4686 absorption also detected in the WHT spectra. 
The detection of He 11 A4686 classifies the primary star in 
HS 1857+5144 as a DAO white dwarf. 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 

2 www.warwick.ac.uk/go/trmarsh 
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Fig. 2. Phase-binned spectra of HS 1857+5144 obtained at McDonald Observatory. These spectra show a clear 
modulation of the emission line strengths with orbital phase. 



2.2. Photometry 

We obtained a total of ~ 60 h of time-series CCD 
photometry of HS 1857+5144 (Fig.Q]) during the period 
July 2003 to June 2006 (Tabled]). Filterless photome- 
try was carried out in 2003/4 using the 1-m Optical 
Ground Station (OGS) and the 0.82-m IAC80 telescope 
at the Observatorio del Teide on Tenerife. Both tele- 
scopes were equipped with Thomson lk x lk pixel CCD 
cameras. Photometric observations were also carried out 
with the 1.2-m telescope at Kryoneri Observatory and 
a 516 x 516 pixel Photometries SI-502 CCD camera. 
The OGS data were reduced in a standard fashion with 
IRAF, and the instrumental magnitudes of the object and 
comparison stars in the field were extracted using the 
point spread function (PSF) packages. Differential mag- 
nitudes of HS 1857+5144 were then computed relative 
to the comparison star 'CI' (USNO-A2.0 1350-10080469: 
R = 13.2, B = 14.7) whose brightness variation was 
found to be negligible against the check star 'C2' (USNO- 
A2.0 1350-10078502: R = 13.5, B = 14.6). The IAC80 and 
Kr yoneri data wer e redu ced using the pipeline described 
by iGiinsicke et al.l (|2004l ) which employs MIDAS for bias 
and dark current subtraction and flat fielding, and per- 
forms apertu re photometry over all visible objects using 
Sextractor ( Bertin fc Arnoutsl 19961 ). Differential magni- 
tudes for HS 1857+5144 were calculated from the Kryoneri 



data using again the comparison star 'CI 'and the check 
star 'C2'. For the IAC80 data, the comparison star 'C3' 
(USNO-A2.0 1350-10078272: R = 14.5, B = 16.1) and 
the check star 'C4' (USNO-A2.0 1350-10079362: R = 16.4, 
B = 17.8) were used. Additional B- and i?-band light 
curves of HS 1857+5144 were obtained at Kryoneri obser- 
vatory in May/ June 2006 and reduced in the same way 
as the filterless data from this telescope. Sample of white 
light, B- and i?-band light curves are displayed in Fig.[H 



Light curve morphology. The light curves of HS 1857+5144 
(Fig.[4j) display a smooth quasi-sinusoidal modulation with 
a period of ~ 6.4 h and peak-to-peak amplitudes of 
0.7 mag in the £?-band, 1.1 mag in the i?-band, and 0.9 mag 
in white light. The minimum in the B-band light curve is 
nearly flat for ~ 0.15 orbital cycle, whereas the shape of 
the minimum in the i?-band is rounder. No sign of the 
typical short-period flickering of accreting systems is de- 
tected, which classifies HS 1857+5144 as a detached bi- 
nary. The low-amplitude scatter seen in the light curves 
in Fig.|4]is caused by residual flat field structures and poor 
tracking of the Kryoneri telescope. The observed periodic 
brightness variation is a characteristic of a large reflec- 
tion effect on the heated face of the secondar y star, irra 
diate d by a ho t primary star (e.g. TW Crv, IChen et al 
19951: KVVel, Iffilditch et all 1 19961: and HS 2333+3927, 
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Fig. 3. Main panel: high resolution spectra of HS 1857+5144 at different orbital phases from 0.02 — 0.77 obtained at 
the WHT. Small windows: close up of the evolution of the H/3 emission line profile in the orbital phase range 0.27 — 0.77 
(from bottom to top). 
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Fig. 4. Filterless, R, and B (from top to bottom) sam- 
ple light curves of HS 1857+5144 obtained with the 1.2-m 
telescope at Kryoneri Observatory. 



3. Analysis 

3.1. Orbital period and ephemeris 

A qualitative inspection of the light curves presented in 
Sect. 12.21 provided an estimate of the orbital period of 
~ 6.4 h. In order to obtain a precise value, we shifted 
the magnitudes of each observing run so that their min- 
ima match that of the 2003 July 21 Kryoneri data, which 
covered an entire orbital cycle. We then subjected the 
combined 2003-2006 dat a to a time-series analysis using 
Schwarzenberg-Czerny's (1996) 0RT method, a variation 
of the analysis-of- variance technique which fits orthogonal 
polynomials to the data folded over a set of trial periods. 
The 0RT periodogram (Fig. [5]) contains an unambiguous 
peak at 3.755 d _1 . A sine fit to the combined photometric 
data defined the following ephemeris: 

T = HJD 2452831.5475(17) + 0.26633357(8) x E , (1) 



Heber et al.l l2004h . Finally, HS 1857+5144 was found at 
a constant mean magnitude of ~ 16.2 throughout our ob- 
serving runs, consistent with USNO-A2.0 measurements 
of HS 1857+5144 (R = 16.3 and B = 15.7). 



where To is defined as the time of inferior conjunc- 
tion of the secondary star (= orbital minimum in the 
light curves). We conclude that the orbital period of 
HS 1857+5144 is T orb = 383.5203 ± 0.0001 min. Figure© 
(bottom two panels) shows the Kryoneri B- and i?-band 
light curves folded according to the above ephemeris. 
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Table 1. Log of the observations. 



Date 


UT 


Telescope 


Filter / 
Grism 


Exp. 

(si 


Frames 


1992 Aug 10 


21:51 


CA2.2 




1500 


1 


2003 Jul 10 


21:03-01:01 


OGS 


clear 


17 


610 


2003 Jul 13 


21:16-01:00 


OGS 


clear 


12 


751 


2003 Jul 21 


19:15-02:25 


KY 


clear 


30 


671 


2003 Jul 23 


18:56-02:32 


KY 


clear 


20-45 


311 


2004 May 21 


01:43-05:25 


IAC80 


clear 


40 


231 


2004 May 23 


00:52-02:07 


KY 


clear 


20-30 


199 


2004 May 25 


21:34-02:26 


KY 


clear 


20 


587 


2004 May 26 


23:14-02:29 


KY 


clear 


20 


335 


2004 May 27 


02:51-05:12 


IAC80 


clear 


40 


286 


2004 Jun 09 


21:40-02:20 


KY 


clear 


20 


608 


2004 Jul 16 


04:13-11:01 


McD 


#43 


600 


31 


2004 Jul 19 


04:25-08:58 


McD 


#43 


600 


4 


2006 May 28 


21:21-02:14 


KY 


R 


70-80 


212 


2006 May 29 


19:53-02:21 


KY 


R 


80 


267 


2006 May 30 


20:58-02:21 


KY 


R 


80 


186 


2006 Jun 04 


20:15-02:14 


KY 


B 


90-110 


196 


2006 Jul 02 


21:17 


WHT 


R1200 


600 


1 


2006 Jul 03 


21:21 


WHT 


R1200 


600 


1 


2006 Jul 04 


21:17 


WHT 


R1200 


600 


1 


2006 Jul 05 


21:16 


WHT 


R1200 


600 


1 



Notes. CA2.2m: 2.2-m telescope, Calar Alto Observatory; 
IAC80: 0.82-m telescope, Observatorio del Teide; KY: 1.2-m 
telescope, Kryoneri Observatory; McD: 2.7-m Harlan J. Smith 
Telescope, McDonald Observatory; OGS: 1-m Optical Ground 
Station, Observatorio del Teide; WHT: 4.2-m William Herschel 
Telescope, Roque de Los Muchachos Observatory. 
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Fig. 5. The DRT periodogram of HS 1857+5144 computed 
from all photometric data. 

3.2. Radial velocities and equivalent widths 

In order to spectroscopically confirm the orbital period of 
HS 1857+5144, we first measured radial velocity variations 
of the H/3, H7, H5, He, and H10 emission lines separately 
by cross-correlating the observed line profiles with a sin- 
gle Gaussian with a full-width at half-maximum (FWHM) 



of 250 km s" 1 for H/3 and He, and of 300 km s" 1 for H7, 
H<5, and H10 within MOLLY. The 0RT periodograms cal- 
culated from the radial velocities of the individual lines 
consistently favoured an orbital frequency of ~ 3.75 d _1 , 
in good agreement with the photometric result. The radial 
velocity amplitudes determined from the different Balmer 
lines varied in the range ~ 175 — 215 km s . 

In order to obtain a more robust measure of the radial 
velocity of the companion star in HS 1857+5144, we de- 
termined an average radial velocity of the Balmer lines by 
fitting simultaneously the Balmer series H/3 to H13 with 
a set of 10 Gaussians. The wavelengths of all Gaussians 
were fixed to their laboratory wavelengths relative to that 
of US, and only the wavelength of US, as well as the 
widths and amplitudes of all 10 Gaussians were used as a 
free parameters. The average Balmer line radial velocities 
are listed in Table[2] and are shown in Fig. [6] (top panel) 
folded over the ephemeris given in Eq.[T] A sine fit to the 
folded velocities and their errors results in an amplitude 
of K em = 185.2 ±4.9 km s" 1 and 7 = -24.0 ± 6.5kms -1 . 
Note that this velocity amplitude does not represent the 
radial velocity amplitude of the centre of mass of the sec- 
ondary star, but that of the illuminated hemisphere. Since 
the centre of light is located closer to the centre of mass 
"of the system than the centre of mass of the secondary 
star, the 'true' radial velocity amplitude of the secondary 
star should therefore be larger than the observed velocity 
amplitude, K cnl = 185.2 ± 4.9kms _1 . We will determine 
a if-correction in Sect. 14. 31 

We also analysed the variation of the equivalent width 
(EW) of the H/3 line (Table©. The Scargle periodogram 
calculated from these measurements contained two equally 
significant signals at 3.409d _1 and 3.749d _1 , the latter of 
which agrees well with the orbital frequency derived from 
the photometry and from the radial velocity variations. 
The equivalent width measurements folded over Eq. |T]) 
are shown in Fig. [6] (second panel from top). As expected 
for an irradiation effect, maximum equivalent width takes 
place at <p ~ 0.5. 

3.3. Photometric and spectroscopic phase relations. 

The assumption that the emission lines in HS 1857+5144 
originate on the inner hemisphere of the secondary star 
as a result of strong irradiation from the hot primary 
star makes specific predictions on the relative phases of 
the photometric and spectroscopic variability. At supe- 
rior conjunction of the secondary star when the irradi- 
ated side faces the observer, tp = 0.5, the system appears 
brightest; the radial velocity of the secondary star is zero 
and crossing from red-shifted to blue-shifted velocities; the 
emission-line strength is at the maximum, and vice versa 
for the inferior conjunction of the secondary star at phase 
zero. Hence, one would expect an agreement in phase 
between the light curve and the equivalent width varia- 
tion, and a 0.25 phase shift betw een those two parameters 
and the radial velocity curve (e.g. Thorstensen et al]|l978 , 



6 



Aungwerojwit et al.: HS 1857+5144: A hot and young pre-cataclysmic variable 



300 
E o 

ffi -200 
10 

2 5 
H 

-5 

-0.5 

0.5 

-0.5 
3 
0.5 



"ft. 



Orbital phase 

Fig. 6. Spectroscopic and photometric data of 
HS 1857+5144 folded over the photometric orbital 
period of 383 min given in Eq.[TJ Top two panels: the 
average of the Balmer radial velocities and H/3 equivalent 
width variations along with the best sine fit (dashed 
line); the filled and open circles represent the Mcdonald 
and WHT data, respectively. The error bars in the radial 
velocity measurements are included in the plot, but of 
similar size as the points. The uncertainties in the values 
of the equivalent widths are dominated by systematic 
effects of the order ~ 1 A. Bottom two panels: i?-band 
and £?-band light curves obtained during May/ June 
2006 along with the synthetic light curves (gray line) 
computed with the program PHOEBE for M wd = 0.72 Mq, 
M sec = 0.21 M , i = 53°, and T wd = 100 000 K. Phase 
zero is defined as inferior conjunction of the secondary 
star. Note that phase of the radial velocity curve is offset 
with respected to equivalent width variation curve and 
light curves by ~ 0.25, consistent with an origin of the 
emission lines on the heated inner hemisphere of the 
secondary star (see Sect. [33] for details). 



1994„ Il996t IVennes & Thorstensenl Il996t lOrosz et all 



199a iHillwig et all 12000" iKawka et all l2002h . Figured 
shows the radial velocity variation of the average of the 
Balmer lines and the equivalent widths of H/3 as well as 
the B- and i?-band light curves folded over Eq.[TJ The 
phase offset between the B- and i?-band light curves with 
respect to the equivalent width variation, as determined 
from sine fits, is ~ 0.004 for the i?-band and ~ 0.018 for 
the B-band. The larger phase offset of the B-band light 
curve is probably related to the fact that it does not cover 
the entire orbital cycle, and hence the sine fit results in 
larger uncertainties. The phase of the folded equivalent 
width variation lags that of the radial velocity curve by 
0.25 ± 0.01 orbital cycle. 



Table 2. The average radial velocities of the Balmer emis- 
sion lines and H/3 equivalent widths of HS 1857+5144 mea- 
sured from the McDonald and WHT spectra . 



HJD 245 


V 


EW 


HJD 245 


V 


EW 




(kms- 1 ) 


(A) 




(kms^ 1 ) 


(A) 


3202.6780 


4.5 + 3.9 


10.8 


3202.8808 


157.4 + 6.7 


4.8 


3202.6857 


-36.2 + 3.4 


9.1 


3202.8884 


161.1 ± 6.2 


5.6 


3202.6954 


-64.0 + 3.6 


9.7 


3202.8954 


161.8 ± 6.4 


6.2 


3202.7037 


-103.4 + 3.5 


8.4 


3202.9058 


140.2 ± 5.5 


7.0 


3202.7134 


-128.7 + 4.1 


8.2 


3202.9134 


115.5 + 5.3 


8.3 


3202.7211 


-150.9 + 5.1 


8.2 


3202.9204 


85.9 + 4.8 


8.1 


3202.7287 


-204.6 ± 6.8 


7.7 


3202.9294 


58.8 + 4.6 


7.9 


3202.7405 


-205.2 ± 6.1 


5.0 


3202.9371 


28.6 + 4.5 


10.0 


3202.7482 


-218.7 + 6.1 


4.4 


3202.9447 


-31.9 + 6.4 


6.8 


3202.7836 




-0.6 


3202.9537 


-44.8 + 4.7 


7.9 


3202.7912 




-2.3 


3202.9614 


-82.5 + 5.0 


9.6 


3202.8013 




-3.0 


3205.6857 


-224.3 ± 8.1 


1.5 


3202.8114 




-3.3 


3205.7468 




-1.9 


3202.8197 




-2.9 


3205.8114 


217.3 ± 10.4 


1.4 


3202.8273 




-3.0 


3205.8760 


28.6 + 4.7 


9.5 


3202.8405 




-0.2 


3919.3887 


-57.8 + 2.8 


8.7 


3202.8482 




-1.0 


3920.3912 


155.6 + 3.5 


4.6 


3202.8572 




0.1 


3921.3888 




-2.7 


3202.8648 


153.1 ± 10.3 


1.2 


3922.3880 


-226.1 +4.7 


2.1 


3202.8725 


140.8 + 8.7 


2.5 









Notes: The statistical error of the H/3 equivalent widths is 
~ 0.01 — 0.05 A, which is negligible. The systematic error, how- 
ever, is of the order ~ 1 A, depending on the details of how the 
continuum flux is determined. 



The phase-dependent behaviour of the emission lines, 
and the relative phases of the photometric, radial velocity, 
and equivalent width variations found in our data corrob- 
orate the hypothesis of the emission lines in HS 1857+5144 
originating on the inner face of the secondary star illumi- 
nated by the hot white dwarf. 

4. Stellar components 

4.1. Light curve solution 

In order to determine additional constraints on the system 
parameters from the observed reflection effect we mod- 
elled the light curves of HS 1857+5144 wi th the 'PHysics 
Of E clipsing BinariEs' program PHOEBeH (jPrsa k, Zwitter 
l2005h. which is built on top of the widely use d WD code 
([Wilson fc Devinnevl[]^7ltlWilsonlll9"79lll990l) . We simul- 
taneously fitted the R- and B-band data obtained at 
Kryoneri Observatory under the following assumptions: 
(a) circular orbits and synchronous rotation of the sec- 
ondary star; (b) stellar surface temperature and bright- 
ness were computed assuming blackbody emission; (c) a 
detailed calculation of the reflection effect was adopted; 
(d) linear limb darkening was chosen, w here the hmbi dark- 
ening coefficient was interpolated from lClaret ( 2000l ): (e) 



http:/ /phoebe.fiz.uni-lj.si/ 
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gravi ty darkenin g exponents of 1 (jvon Zeipell [1924) and 
0.32 (|Lucvlll967l) were used for radiative and convective 
stars, respectively; (f) no contribution fluxes from a spot 
or third light were applied. 

In our analysis we tested a wide range of white dwarf 
mass, covering M wd = 0.3 — 1.4M . We then assumed 
an M-type companion star, testing the whole range of 
spectral type M9-M0V, corresponding to masses, radii, 
and temperatures of M scc ~ 0.07 — 0.53 Mq, R scc ~ 
0.11 - 0.56 R Q , and T scc ~ 2300 - 3800 K, interpolated 
from Rebassa-Mansergas et al. (in preparation). Earlier 
spectral types than M0 V would imply extremely massive 
white d warf s ( S ect . 1 4 . 3 P which are excluded by the spectral 
fit (Sect. l4~2"|i 4 l . This approach allows us to search for pos- 
sible solutions over a large range of possible mass ratios, 
q = M sec /M w( j — 0.05 — 1.77. For each input q we fixed 
-Porb = 0.26633 d and T scc according to the selected spec- 
tral type of the companion star. The following parameters 
were free in the fits: q, white dwarf temperature (T wd ), 
orbital inclination (i), surface potentials for both compo- 
nents, and albedo of the secondary star (ALB2). The fits 
for early type companions, M2-M0V, do not reproduce 
well the observed B and R light curves, which supports the 
exclusion of early-type donors outlined above. We found 
that, independently from the details of a given fit, the sys- 
tem must contain a hot white dwarf with T wd > 60 000 K 
and a cool component with an albedo higher than that 
of a normal M star (ALB2 > 0.5) to reproduce the large 
amplitude observed in the light curves. Such high T w d is 
also confirmed by the spectral fit to the WHT faint-phase 
spectrum in the following section. Fitting the B- and R- 
band light curves alone provides a fairly large range of 
possible system parameters, M w d — 0.3 — 1.4M Q , M scc ~ 
O.O66-O.367M , (spectral type M9-M3V), i ~ 40°-60°, 
and T wd - 60 000 - 100 000 K. For a given M soc , a more 
massive white dwarf requires a larger inclination and a 
higher T wd . 

Figure[6] (bottom two panels) presents the correspond- 
ing R- and i?-band synthetic light curves from the pro- 
gram PHOEBE for M wd = 0.72 M , M aec = 0.21 M Q , 
i = 53°, and T wd = 100 000K, along with the observed 
light curves folded over the ephemeris in Eq.[TJ The choice 
of this particular set of parameter is detailed below in 
Sect. 14.41 but other fits in the parameter range given above 
fit the data equally well. 



1.2 




H(3 



0.8 



-100 





AA [ A] 



Fig. 7. H/3 and H7 from the WHT faint-phase spectrum 
((f = 0.02) fitted with an LTE model for T wd = 100 000K 
and log g = 8.0. 



of M w d and T w d, using both a grid of LTE pure-hydrogen 
models ([Koester et al.1120051) . and a grid of NLTE models 
with a variety of He abundance^]. We fitted the H/3 and 
H7 absorption lines after normalising the continua of the 
observed data and the model spectra in the same way us- 
ing a third-order polynomial. The fits suggest 80 000 K ^ 



T wd < 100 000K and 7.5 <, log g <, 8.5 (corresponding 
to 0.6 M & M wd l.QM Q ) for the LTE models, and 
70 000 K < T wd < 100 000 K and 8.0 ^ log g < 8.5 (cor- 
responding to 0.7M Q < M wd < l.OM ) for the NLTE 
models. Figure[7] shows the best LTE fit to H/3 and H7 for 
T wd = 100 000 K and log g = 8.0. These numbers should 
be considered as rough estimates only, as the optical spec- 
trum of HS 1857+5144 is contaminated by flux from the 
companion star, which is very difficult to quantify. A more 
reliable temperature and mass estimate would require far- 
ultraviolet data, where the white dwarf dominates the 
emission of the system ([Good et al.l 120041) . We conclude 



from this qualitative spectral analysis that the white dwarf 
in HS 1857+5144 is indeed very hot and, given the high 
surface gravity preferred by the fits, is more likely to be a 
white dwarf than a subdwarf. The detection of He 11 A4686 
then qualifies the primary as a DAO white dwarf. 



4.2. Spectral fit 

We performed a spectral fit to the faint-phase WHT spec- 
trum of HS 1857+5144 to obtain an independent estimate 



4 A strict upper limit on the mass of the companion comes 
from the fact t hat it is not R oche-lobe filling. Using p ~ 
107P o ^(h) fe.g. lEggletonlllOS:! ), with p, the average density 
of the donor, and P or b(h), the orbital period in hours for a 
Roche-lobe filling star, the maximum mass of a main-sequence 
companion in HS 1857+5144 is M Bec < 0.72 M Q , corresponding 
to a spectral type K3 V or later. 



4.3. K-correction and mass ratio-inclination 
constraints 

As mentioned already in Sect. 13.21 the emission lines in 
HS 1857+5144 trace the orbit of the centre of light of the 
illuminated hemisphere of the secondary star, and not its 
centre of mass. Hence, for a dynamic assessment of the 
binary parameters, the measured velocity K em (Sect. l3~2"|) 
has to be corrected accordingly. The radial velocity ampli- 



5 http:/ /astro. uni-tuebingen.de/~rauch 
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tude of the secondary sta r's centre of ma s s, K* prrr , r can 
be expressed according to Wade fc Home! ( 19881 ) as 



a;, 



K, 



l-(l + q)(AR/a) 



where AR is the displacement of the centre of light from 
the centre of mass of the secondary star, with ^ AR < 
R scc (AR = implies that the centre of light coincides 
with the centre of mass of the secondary star, whereas 
AR = R scc is the maximum possible displacement, where 
all the light comes from a small region on the secondary 
star closest to the primary star). Assuming that the emis- 
sion due to irradiation is distributed uniformly over the in- 
ner hemisphere of the secondary sta r, and zero on its unir- 
radiated face, AR = (4/3ir)R aar - (Wade fc Horn! Il988t 
Wood et al.lll995HOrosz et alJll999t IVennes et alJll999l ). 
The expected radial velocity of the secondary star 

(-A seC5Ca i) is 



K. 



scc,cal 



27ra sini 

Porb(l + g) 



where a is the binary separation. Equating A" sccxor = 
Agcccai then gives a unique q for a given choice of i, 
Hence, a fixed value of i projects onto a one-dimensional 
curve within the (M w( j,M sec ) plane, and for the possi- 
ble range of parameters considered here, those curves are 
nearly straight lines. 

4.4. Combined constraints 



In Sect. 14. 1114. 31 above we have outlined what type of con- 
straints on the system parameters of HS 1857+5144 can 
be derived from the observed light curves, radial velocity 
variations, and the spectrum of the primary star. Here, we 
will combine all those independent constraints. 

In a first step, we impose a range of i = 40° — 60°, as 
suggested by the set of PHOEBE fits to the B- and i?-band 
light curves, on the combinations of (M w d,M sec ) which 
are consistent with the corrected radial velocity -ftT SO c,cor of 
the secondary star. The resulting parameter range is indi- 
cated by the gray shaded area in Fig. [8] In a second step, 
we inspected the individual light curve fits from the grid 
of PHOEBE runs, and require the inclination of a model for 
a given (M w d,M scc ) to fall within ±5° of the correspond- 
ing inclination constraint from the radial velocity of the 
secondary star. We introduce this "fuzziness" in inclina- 
tion as a measure to account for systematic uncertainties 
within the if-correction and the light curve fits. Possible 
combinations of (M w( j,M sec ) are indicated by filled cir- 
cles in Fig.[8j and trace a somewhat narrower band than 
the initial i = 40° — 60° constraint. A final constraint 
comes from the spectral fit of the WHT faint phase spec- 
trum (Sect.O, which implied 0.6 M <, M wd <, 1.0 M , 
shown as vertical dashed lines in Fig.[S] 

The combination of all constraints suggests 0.15 M 
M scc < 0.30 M Q (spectral type M6-M4V), 0.6 M ^ 
M wd < 1.0 M©, 70 000 K < T wd < 100 000 K, and 45° SS 




0.4 0.6 



M WD [MJ 

Fig. 8. Photometric and spectroscopic constraints on 
(M wd , Ms ec ). Gray shade area represents possible dynam- 
ical solutions from the 'K-corrected' radial velocity of the 
secondary star for K em = 185 kms^ 1 and for i = 40°— 60°. 
Filled circles represent possible solutions from the light 
curve analysis for given (Af wd , M SQC ) in which i agrees with 
spectroscopic constraint within ±5°. Dashed lines are up- 
per and lower limits on M wc \ derived from the spectral fit 
(see Sect. [O] for details). 
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Fig. 9. The BRJHK S fluxes of HS 1857+5144 (filled tri- 
angles) at orbital maximum (ip ~ 0.5) along with an 
example of a plausible fit (solid line) with the sum of 
a white dwarf spectrum (dots) and a blackbody (dash 
line) representing the heated side of the low-mass com- 

70 000 K, i? wrl = 1.3 x 10 9 cm, 



panion, assuming T w( j 
1 cm, 

460 pc (see Sect.[0]for details) 



at a distance of 



i 55° . A substantial improvement on this set of param- 
eters will require measuring the radial velocity amplitude 
of the white dwarf, K^d, and determining an accurate 
temperature for the primary. Both types of measurements 
could be easily obtained from time-resolved ultraviolet 
spectroscopy. 
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4.5. 2MASS magnitudes 



HS 18 57+5144 is detected in the 2MASS (jSkrutskie et al 
20061) at J = 15.09+0.04, #15.05+0.08, K s = 14.76+0.14. 
Our ephemeris is good enough to establish the orbital 
phase of the 2MASS data, which is very close to or- 
bital maximum, ip ~ 0.5. In order to test to what extent 
the optical-infrared spectral energy distribution (SED) 
of HS 1857+5144 is compatible with the stellar photom- 
etry and spectroscopy, we modelled the BRJHK S pa- 
rameter range obtained from the analysis of the time- 
resolved ma gnitudes with the su m of a white dwarf spec- 
trum from ( Koester et al. 2005f ) and a blackbody repre- 
senting the contribution of a (heated) low-mass compan- 
ion. The data leaves some freedom in the exact parame- 
ters, however, as an example, a fairly good fit is achieved 
for T w d = 70 000 K, i? w d = 1.3 x 10 9 cm, corresponding 
to M w( j ~ 0.6 M Q , i? sec = 1.3 x 10 10 cm, corresponding to 
spectral type of M6, and T sec = 6000 K, at a distance of 
460 pc (Fig. [9]). This solution coincides well with the pho- 
tometric and spectroscopic constraints in Fig.[H Overall 
our simple model confirms the stellar parameters estab- 
lished in the previous sections, with a slight preference for 
a white dwarf mass in the range Af w d = 0.6 — 0.8 M , 
and a companion star with a radius near the lower end 
of the determined range, i.e. corresponding to a spectral 
type M6-M5, and a distance to the system of 290-460 pc. 



5. Discussion 

The analysis presented in Sect. 14. 21 suggests that 
HS 1857+5144 contains a hot white dwarf with T w d — 
70 000 - 100 000K. T he implied cooling ag e of the white 
dwarf is 1.2— 6xl0 5 yr ( Bergeron et al.ll995 . and Bergeron 
2002, private communication), making HS 1857+5144 one 
of the youn gest PCEBs known so f ar. F ollowing the pre- 
scription of Schreiber fc Gansicke ( 20031 ). and assuming 
the range of system parameters established in Sect. [4] as 
well as "classical" magnetic braking for the angular mo- 
mentum loss mechanism, we estimate the period at which 
HS 1857+5144 left the common envelope phase to be 
Pce — 0.266334 — 0.266345 d, very close to its present 
orbital period. HS 1857+5144 will evolve within the next 
~ 0.4 — 1.3 x 10 10 yr into a semidetached CV configu- 
ration, and start mass transfer at an orbital period of 
~ 0.08 — 0.13 d, i.e. within or below the period gap. The 
large uncertainties on the future evolution are a conse- 
quence of the limited constraints on the system param- 
eters. Additional systematic uncertainties in the actual 
strength of angular momentum loss from the orbit have 
not been taken into account. 

Among ~ 40 prev i ously know n PCEBs liste d 



m 



Schreiber fc Gansickd (|2003l ). and IShimanskv et al 



( 20061 ). only 8 systems display a reflection effect compara- 



probes into our understanding of common envelope evo- 
lution. 

A large reflection effect is expected for those PCEBs 
containing a hot subdwarf, because of the larger luminos- 
ity compared to a white dwarf of the same temperature, 
and indeed the majority of known PCEBs with a large 
reflection effect have sdO primary stars. So far only one 
PCEB, HS 1136+664 6, containing a ho t white dwarf simi- 
lar to HS 1857+5144 (jSing et al.ll2004h is known. The sec- 
ondary star in HS 1136+6646 has been suggested to be 
a K7-4V star on the basis of its spectral type, which 
appears too early for the estimated mass of 0.34 M Q . 
However, ISing et all (|2004h discuss the possibility that the 
secondary is overluminous as it is still out of thermal equi- 
librium after accreting at a high rate during the common 
envelope phase. The amplitude of the reflection effect in 
HS 1136+6646 is much lower than in HS 1857+5144, which 
is consistent with its longer orbital period of 0.84 d. The 
other system most similar to HS 1857+5411 is BEUMa, 
which has been classified as a borderline o bject between 



an sdO subdwarf and a DAO white dwarf (Lieb ert et al 



ble to that of HS 1857+5144 (Tabled . All those systems 
contain extremely hot white dwarfs or subdwarfs, and all 
are very young PCEBs that may serve as observational 



1995 : Ferguson et al. 1999), and is associated with a plan- 
etary nebula. It is interesting to note that six out of the 
nine systems listed in Table[3]are still embedded in a plan- 
etary nebula. Our long-slit spectroscopy of HS 1857+5144 
does not reveal any noticeable trace of extended emission 
around Ha, though a deep Ha image testing for remnant 
nebular emission would be useful. Similarly, no sign of ex- 
tended Ha emission around HS 1136+6646 has been ob- 
served. While the majority of very young PCEBs are still 
embedded in their planetary nebulae/common envelopes, 
the discovery of HS 1857+5144 and HS 1136+6646 sug- 
gests that the envelope can be dispersed rather quickly. 

6. Conclusions 

We have identified a pre-CV, HS 1857+5144, containing 
a hot DAO white dwarf and a low mass MV star with 
an orbital period of P rb = 383.52 min. The light curves 
of HS 1857+5144 display a very large reflection effect with 
peak-to-peak amplitude of 0.7 and 1.1 mag in the B and R 
bands, respectively. Combining the results of our spectro- 
scopic and photometric analysis, we constrain the system 
parameter to 0.15 M Q ^ M scc Si 0.30 M Q (spectral type 
M6-M4V), 0.6 M ^ M wd <, 1.0 M©, 70 000 K < T wd ^ 
100 000 K, and 45° & i & 55°. Taking these parameters at 
face value HS 1857+5144 is one of the youngest PCEBs 
known so far and has just emerged from the post common 
envelope phase. The system will eventually evolve into a 
cataclysmic variable, and start mass transfer within or be- 
low the 2 — 3 h orbital period gap. 
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Table 3. PCEBs with a large reflection effect. 



Object 


Porb 


SP1 


SPsec 


Ti 


Mi 




Reflection effect [mag] 


Ref. 




[d] 






[K] 


[M ] 


[M ] 


B 


V 


R 




NNSer 


0.130 


DA 


M4.75V 


57 000 + 3000 


0.54 ± 0.05 


0.150 + 0.008 


0.33 


0.49 


0.772 


1, 2 


HS 1857+5144 


0.266 


DAO 


~M6-M4V 


~ 70 000 - 100 000 


~ 0.6 - 1.0 


~ 0.15 -0.30 


0.7 




1.1 


3 


TWCrv 


0.328 


sdO 


MV 


105 000 ± 20 000 


~ 0.55-0.61 


< 0.3 


0.74 


0.85 


0.93 


4, 5 


KVVel 


0.357 


sdO,PN 


MV 


77 000 + 3000 


0.63 + 0.03 


0.23 + 0.01 


0.49 


0.55 


0.61 


6 


V477Lyr 


0.472 


sdOB,PN 




60 000 + 10 000 


0.51 + 0.07 


0.15 + 0.02 


0.5 


0.6 




7, 8 


V664 Cas 


0.582 


sdO,PN 


K5-F5V 


83 000 ± 6000 








1.15 




9, 5 


VWPyx 


0.676 


sdO,PN 




85 000 ± 6000 








1.36 




10, 11 


Abell65 


1 


sd?,PN 




~ 80 000 








> 0.5 




12, 13 


BEUMa 


2.291 


sdO/DAO,PN 


K4-3V 


105 000 ± 5000 


0.70 + 0.07 


0.36 + 0.07 




~ 1.3 




14, 15, 16 



References: (1 ) iHaefnerl (|l989h 



(2) lHaefner et af 1 



20041): (3) this w ork; (4) IChen et all (|l995l ): (5) lExter et all d2005l): 



iHilditch etahl dl996l); ( 7) bo nd+grauer-87;_ ( 8)|PoUacco fe Bell dl994l); (9) shima nskiie tal04-l; (10) iKohoutek fe Schnur! (1198: 



(11) lExter et all d2003h: (12 ) iBond fe Livid l| 19901 ): (13) I Walsh fc Waltonl (|l996h : (14) iFerguson et all (|l987l ): (15) IWood etahl 
|l995h : (16) lFerguson et all l| 19991 ) 



4f } 
82); 
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